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In order to clarify common assumptions on the form of energy and momentum in elas- 
ticity, a generalized conservation format is proposed for finite elasticity, in which total 
energy and momentum are not specified a priori. Velocity, stress, and total energy are 
assumed to depend constitutively on deformation gradient and momentum in a manner 
restricted by a dissipation principle and certain mild invariance requirements. Under 
these assumptions, representations are obtained for energy and momentum, demonstrat- 
ing that (i) the total energy splits into separate internal and kinetic contributions, and 
(ii) the momentum is linear in the velocity. It is further shown that, if the stress re- 
sponse is strongly elliptic, the classical specifications for kinetic energy and momentum 
are sufficient to give elasticity the standard format of a quasilinear hyperbolic system. 

Keywords: finite elasticity, hyperbolic systems of conservation laws. 

1. Introduction 

In continuum mechanics the total energy of a body part is generally assumed to be 
the sum of an internal energy plus a kinetic energy, with the latter quadratic in the 
velocity; moreover, the momentum is assumed to be linear in the velc 
expressions for total energy and momentum are, with few exceptions,! 
accepted without further discussion. One may ask, however, whether they could be 
derived from first principles. Such a derivation would perhaps allow for alternative 
expressions, and yet lead to a plausible mechanics; or, while only conceptually im- 
portant for well-established continuum theories, it would suggest a modus operandi 
in the case of continuum theories where identification and representation of inertial 
quantities is by no means obvious {e.g., theories of liquid crystals, deformable fer- 
romagnets, moving phase boundaries, crack dynamics, adhesion and peeling, . . . ). 
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Recently, a derivation based upon invariance of the internal power and linearity of 
the inertial power has been proposed for classical Cauchy continuajL3 the format of 
this derivation is sufficiently robust to allow for various generalizations.ErET E3iiJll3 
Here we follow an alternative course: we base our approach on a presumption of 
mathematical structure, namely, that the representations for energy and momen- 
tum in finite elasticity be such that the resulting evolution problem have the form 
of a weakly- hyperbolic system of conservation laws, compatible with a dissipation 
inequality and certain mild invariance requirements. 

In Section 2 we first recall briefly the formulation of initial-value problems for 
general quasilinear hyperbolic systems of conservation laws.B'HLj Within this formu- 
lation, the standard format for nonlinear elasticity assigns the role of state variables 
to the deformation gradient F = d x f and the velocity v = dtf in a motion f(x, t); 
further, the momentum is taken proportional to velocity, the stress S = d-pa is 
determined by the stored energy <x(F), a and the total energy r is stored energy plus 
kinetic energy, with the latter proportional to |v| 2 . 

We introduce next a generalization of this conservation format for elasticity, in 
which the total energy, the stress, and the velocity are given by general constitutive 
prescriptions r = f(F,p) etc. in terms of a pair of kinematic state variables, the 
tensor F and the vector p, to be interpreted as deformation gradient and momentum 
but initially unrelated to motion. As a criterium of constitutive admissibility for 
f(F,p), S(F,p), and v(F, p), we postulate that the purely mechanical dissipation 
inequality 

f(F,p) < S(F,p).F + v(F,p).p 

be identically satisfied in all admissible smooth processes consistent with the con- 
servation laws 

F(ac,t)-Div(v(x,t)<8)l) = 0, p(x,t) - Div S(x,t) = 0. (1.1) 

In this formulation of evolutionary elasticity, momentum replaces velocity as a kine- 
matic state variable. No relation of momentum to mass and velocity is postulated; 
in particular, there is no need to assume that the referential mass density be con- 
stant in order to apply the general conservation format to elasticity. The format 
we propose requires that quantities describing boundary fluxes be the subject of 
constitutive prescriptions; we then introduce such prescriptions for both stress (= 
boundary flux of momentum) and velocity (= boundary flux of deformation gradi- 
ent). 

The purpose of our generalized formulation of elasticity is to determine the 
constitutive specifications for energy, momentum, and stress that are compatible 
with both the dissipation inequality and the conservation laws. Section 3 provides 
these results. To obtain them, we require arbitrary local continuations (F,p) of 
a kinematic process through (F,p); since we do not include source terms in the 
conservation laws, we can only count on initial conditions. Two technical conditions, 



a In elasticity, a purely mechanical context, stored energy replaces internal energy. 
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together akin to hyperbolicity but weaker, are to be satisfied by the constitutive 
mappings for velocity and stress. The first condition is 

• normality of v(F,p) with respect to p for each fixed F. 

A normal v(F,p) has a local inverse p(F,v) for each fixed F; it follows that the 
stress mapping S(F, v) = S(F,p(F,v)) is well-defined. The second condition is 

• ellipticity of S(F,v) with respect to F for each fixed v. 

Under these assumptions, arbitrariness in local continuation is guaranteed, and 
thermodynamically admissible constitutive mappings are shown to satisfy 

v(F,p)=9 p f(F,p), S(F,p)=G>Ff(F,p). (1.2) 

These relations are necessary precursors to the representations for velocity and 
energy, which are obtained under two additional constitutive assumptions: 

• Galilean variance of velocity in a translational change in observer, 

• parity of total energy with respect to momentum. 
These lead to the representations: 

v(F,p)=Vp, f(F,p) = fc(p) + *(F), £(p) = ip.Vp, (1.3) 

with V a symmetric, invertible, time-independent tensor. We conclude that, for 
elasticity to have the mathematical structure of a system of conservation laws com- 
patible with a dissipation inequality, momentum and kinetic energy must have the 
representations 

p = Mv, «=iv-Mv, (1.4) 

where the symmetric mass-density tensor M is the inverse of V; moreover, since M 
inherits time-independence from V, mass must be conserved. These representations 
are similar to the classical ones, b but slightly more general due to the tensorial 
nature of mass density. 

Relations (1.3)2 an d (1.4) — the main result of this paper — are the same arrived 
at when the evolutionary problem for general Cauchy continua is put in the invari- 
ance format introduced in Ref. 13. Remarkably, the assumptions implying (1.4) are 
compatible with hyperbolicity of the system (1.1) of conservation laws if the mass- 
density tensor has the classical form M = pi (with p > 0) and the stress mapping 
is strongly elliptic. 

'-'Namely, the representations 

p = p v and k = — p v ■ v , 



with p the referential mass density. Classically, one also tacitly assumes the additive splitting (1.3)2 
of t hj^ totaJ e nerg-v as well as that the stress is the derivative of the stored energy, a consequence 



of ([ljh and (|L3|)2 
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2. Elasticity in Conservation Format 

2.1. General hyperbolic systems with involutions 

The general format of theinitial- value problems for quasilinear hyperbolic systems 
of conservation laws isBBO 

m 

d t U{x,t) + J2 d <*G a (U{x,t)) = 0, (i,t)el m x(0,oo), (2.1) 

a=l 

U{x,0) = U {x), xgl m . (2.2) 

Here dt stands for d/dt and d a for d/dx a ; U, the state vector, takes values in a 
subset O of W 1 . The m constitutive maps G a from O into R n are supposed to be 
smooth and such as to satisfy the following hyperbolicity condition: 

(H) for each fixed U and for each unit vector w in R m , the (n x n)-matrix 

Y, w -d v G a (2.3) 

a=l 

has real proper numbers and n linearly independent proper vectors. 

Generally, the systems of conservation laws of interest for applications have an 
associated "entropy" -"entropy flux" pair (77, q), with the smooth mappings r](U) 
and q(J7) defined over O and taking values in K and R m , respectively. The role of 



this pair of mappings is to characterize the admissible solutions to ( |2.l|) -(2.2) as 
those satisfying the generalized "entropy condition" 

m 

d t n(U) + J2daq a (U) < 0. (2.4) 

a=l 

Often the conservation system ( J2.1J ) has additional geometrical structure em- 
bodied in a system of k "involutions, "u that is, of k linear differential equations 

m 

A a d a U = 0, A a — a constant (k x n)-matrix, (2-5) 



to be satisfied by the solutions of (2T)-(2 : S) so long as they are satisfied by the 



initial data (2^). In the presence of involutions the hyperbolicity condition generally 
needs to be reformulated; in particular, the number of linearly independent proper 
vectors may be less than n. 

2.2. Finite elasticity: the standard format 

The system of finite elasticity in the referential formulation is a standard example 
of a hyperbolic system with involutions (provided the mass density is constant). To 
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show that t he general format just described applies, consider the elasticity system 



written as 



F — Vv = 0, {compatibility condition) (2-6) 
v — Div S = 0, [momentum balance) (2-7) 

where, for convenience, the referential mass density is presumed to have unit value 



(whence the coincidence of momentum with velocity in (2/7)), and where S is in- 
terpreted as the Piola stress tensor. c 

The state vector U is the pair (F, v) of a tensor F and a vector v. To insure that 
these can be interpreted, respectively, as the deformation gradient and the velocity 
in a motion, that is, to insure that there is a motion f(x,t) such that F, v, and / 
satisfy 

F(x, t) = Vf{x, t) and v[x, t) = d t f(x, t), (2.8) 



it is required that F obey the involution following from (2.8)i, namely 



(V(Fa))b = (V(Fb))a (2.9) 

for all constant vectors a, b (cf. ([2~5|)). d 

As to the "entropy" - "entropy flux" pair, one chooses 

7y(F,v) = «(v)+«7(F), (2.10) 

q(F,v) = -S T (F)v, (2.11) 

with 

«(v) = ivv, (2.12) 

the kinetic energy, er(F) the stored energy, and 

S(F) =9 F cr(F). (2.13) 

Thus the "entropy" in this case is the total energy (re + a) per unit volume, the 
"entropy flux" is the flux necessary to balance the total energy, and the generalized 



"entropy condition" (2.4) takes the simple form 

d t (n + o-) -Div (S T v) = 0. (2.14) 

Remarkably, this relation is consistent with the Second Law of thermodynamics, 
when the latter is formulated in the form of the Clausius-Duhem inequality. 

c As is common in continuum physics, we have used a superposed dot in place of dt and the nabla 
symbol in place of d x . 

d As is w ell-known, given the pair (F,v), the kinematical compatibility conditio ns ( p.fij ) 



and ( p.9[ ) are sufficient for the existence of a local solution / of the system (2.8). 
Here and henceforth, unless otherwise specified a vector is meant to be an element of a three- 
dimensional inner-product vector space V, and a tensor is a linear transformation of V into itself. 
Moreover, whenever needed or simply convenient, V will be thought of as being endowed with an 
orthonormal basis {c a , a = 1,2,3}. All components of vectors, tensors, etc. will be tacitly taken 
with respect to this basis and related constructs. 
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2.3. Finite elasticity: a generalized format 

Consider now the following generalization of the standard conservation format for 
finite elasticity. 

For the state vector U = (F, p), and for 1 the identity tensor, let the initial- value 



problem (2.1)— (2.2) take the form 

F(x, t) - Div (v(x, t) <g> 1) = 0, (2.15) 

p(x,t) -DivS(a:,t) = 0, (2.16) 

F(x,0) = F„(a:), p(x,0) = Po (x). (2.17) 



Next, as the generalized format suggests, let the conservation laws ( 2.1 5| ) and ( 2.1 6| ) 



for the state variables "deformation gradient" F and "momentum" p be supple- 
mented by constitutive prescriptions for velocity and stress as functions of state: 

v = v(F,p), (2.18) 
S = S(F,p). (2.19) 

Finally, let the generality of these constitutive prescriptions and of the additional 
prescription for the total energy: 

r = f(F,p), (2.20) 

be restricted by requiring consistency with the Second Law, under form of the 
dissipation inequality 

f (F, p) < S(F, p) • F + v(F, p) • p, (2.21) 

to be identically satisfied in all admissible smooth processes. 

As to our present formulation of finite elasticity, some comments are in or- 



der. Firstly, the state variables for system (2.15)— (2.16) are a tensor F and a 
vector p. Formally, neither F need be the gradient of a deformation (this geometric 
compatibility problem may be solved once a solution of the initial-value problem 



(2.15)— (2.17) has been found, but at this stage we do not require that F satisfy the 
involution ( P-9| )) nor that p need have any a priori specified relation with mass den- 
sity and velocity. The compatibility condition (|2.6|) is generalized to the evolution 



equation ( |2.15 ), which we interpret as the balance of F by way of the boundary- flux 



tensor v ® 1. Similarly, the momentum balance (2.7) is generalized to (2.16), the 
balance of p by way of the boundary-flux tensor S. Less formally, we regard (2.16) 
as a force balance, with — p the inertial force (the only external force considered 
here) and Div S the internal force per unit referential volume. 

Secondly, just as the momentum is not given any a priori representation, the 
other inertial manifestation — kinetic energy — is not given any a priori representa- 
tion as well. In fact, neither the total energy is split into a kinetic part and an 
internal part nor the Piola stress is introduced by differentiation of the latter. 

e In all constitutive mappings we leave tacit a possible explicit dependence on the space variable, 
but exclude any explicit time dependence. 
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Our next step is to introduce and exploit a set of reasonable hypotheses allowing 
to conclude that slight generalizations of the classical constitutive specifications for 
energy, momentum, and stress, are the only specifications compatible with both the 
conservation laws ( p. 15 ), (2.16) and the dissipation inequality (2.21). 



3. Representations of Energy and Momentum 

3.1. Normality and ellipticity. Thermodynamical admissibility 

As is customary in constitutive theories of continuous media after a famous paper by 
Coleman & Nolljd the dissipation inequality restricts the choice of admissible repre- 
sentations of the constitutive mappings. These restrictions are derivable under the 
crucial assumption that any accessible kinematic state (F, p) can be reached by a 
path that has an arbitrary local continuation (F,p). Such indispensable arbitrari- 
ness is usually taken to be granted by the presence in the basic balance equations 
of source terms regarded as controls at our disposal. Since sources are here set to 
null, we revert to the other piece of data that in principle can be made to vary 
arbitrarily, that is, initial conditions. f Our first proposition establishes that such 
an indispensable arbitrariness is indeed granted under reasonable conditions on the 
balance equations and the constitutive mappings for velocity and stress. 

The constitutive mappings v and S are taken to be continuously differentiable 
and to satisfy two conditions, one each. The first condition applies to v: 

(N) (Normality) For each F fixed, the map p ^ v(F, p) is normal, i.e., the matrix 
N, defined by 

N ih (F,p) = d Ph Vi(F,p), 

is invertible for all p. 

This condition formalizes the expectation that momentum and velocity be locally 
in one-to-one correspondence, in the sense that the map (F,v) i— » (F, p = p(F,v)) 
is well-defined, with 

v(F,p(F,v)) =v. 

Let now 

S(F,v):=S(F, P (F,v)). (3.1) 
The second condition, rather than directly to S(F, p), applies to S(F,v): 

(E) (Ellipticity) For each v fixed, the map F i— > S(F,v) is elliptic in the sense of 
Petrovsky, i.e., the matrix E, defined by 

Eih{F,v,a) — dF hk Sij(F,v) ajdk for each vector a, 

is invertible for all F. 

This idea is not new BEET although we know neither of any previous use in the context of 
elasticity nor of anything more than a generic claim of local existence for the initial- value problem 
under examination. 
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Normality of v(F, •) and ellipticity of S(-,v) do not imply hyperbolicity of the 
system (2.15)-( p.l6 ). As discussed later in this section, some strengthening of these 
assumptions is needed to obtain hyberbolicity. For the time being, we consider the 
following initial-value problem: 



For each T > and if! 3 fixed, find 1 i— > (F(x, t), p(x, t)) such that 
t(x,t) = Vv(x,t), 



p(x,t) 

for t 6 (0, T) and, moreover, 

F(x,0) 
v(s,0) 



DivS(F(x,t),v(x,t)), 

A + a-(x - x D )(b ® a), 
B(x - x ) + c. 



(3.2) 
(3.3) 

(3.4) 
(3.5) 



Here x D is a given point of R 3 , A and B are two given matrices, and a, b, and 



c three given vectors. This problem is obtained by supplementing (2.15)— (2.16) by 
the constitutive relation (|3.l[) and the initial conditions (3.4)-(p^), which involve 
the state variables (F, v); note that, due to (N), (|3.4|)-(3.5) induce a unique pair of 
initial conditions for the state variables (F, p). 

Proposition 1 (Arbitrariness of initial time-rates of the state variables) 

Assume that the constitutive mappings v and S are, respectively, normal and ellip- 
tic. Moreover, assume that, for each fixed triplet (A, a, c), the initial- value problem 



(3.2)— (3.5) has a classical solution up to time T for all choices of B, b, and x . Then 
the initial time-rate (F(x , 0), p(x D , 0)) at x a of the state vector can be assigned 
arbitrary values. 



Proof. Note firstly that the evolution equations (|3.2|)-(3.3) can be written as 



F-ij Oj Vi , 

t>i = 9 Fhk Sij (F, v) djFhk + d Vk Sij (F, v) djv k . 

Moreover, the initial conditions ( |3.4[ )-( |3~5| ) yield 

(F(x o ,0),v(x o ,0)) = (A,c), 
(djFhk(x o ,0),9jV k (x o ,0)) = (b h a k aj,B kj ). 



(3.6) 
(3.7) 

(3.8) 
(3.9) 



Then, since the initial-value problem ( |3.2j )-(3.5) has classical solutions, with the 
notations introduced in the statement of condition (E) we have that 



F%j (x , 0) 
Pi(x o ,0) 



Bij , 

d Vk Sij(A, c) B kj + E ih (A, c; a) b h . 



(3.10) 
(3.11) 



When B and b are arbitrarily varied, the invertibility of E(A, c; a) guarantees that 
the second term on the right of the last relation and, hence, p(x , 0) take arbitrary 
values, not only F(x o ,0). □ 
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Proposition 2 (Thermodynamically admisssible constitutive relations) 

Under the assumptions of Proposition 1, necessary and sufficient conditions for the 



dissipation inequality (2.21) to be satisfied in all admissible smooth processes of the 



materials described by the constitutive relations (2.18)— (2.19) are 

v(F,p) = 9 p f(F,p), (3.12) 
S(F,p) = a F f(F,p). (3.13) 

Moreover, if f is twice continuosly differentiable, then 

d p S = 5 F v. (3.14) 



Proof. From (2.21) and by the chain rule, 

(9pf-S)-F+(a p f-v)-p < 0. (3.15) 

Sufficiency of (3.11)— (3.12) follows by substitution. Since (F,p) can be chosen 
arbitrarily, at least initially (Proposition 1), we obtain necessity. □ 

3.2. Galilean variance and parity. Representation results 

It is common in mechanics to impose some type of invariance requirements on the 
constitutive relations. Here we consider a notion of translational change in observer 
appropriate to our present choice of the state vector. Since neither F nor p are 
presumed to have any relation to the motion, we cannot deduce their variance laws 
from the standard notion of observer change. Instead, we say that the pairs (F,p) 
and (F*,p*) are related by a translational change in observer whenever there is a 
vector d, independent of time, such that 

F^F* = F, 

p^p' = p + d. (3.16) 

Thus, while F is invariant, p is not; d measures the "invariance defect" of the latter. 
Note that when the momentum p is given the classical representation p = pv, this 
notion reduces to the standard one, provided one chooses d = pe, with e a constant 
vector. As to the behavior in a translational observer change of the constitutive 
mapping delivering the velocity, we introduce the following assumption. 

(G) (Galilean Variance of Velocity) In a translational change in observer, 

v(F, p + d) v(F, p) = c v (d) for each vector d, (3.17) 

where the mapping c v specifies the invariance defect of the velocity mapping 
v. 

Our next assumption specifies the parity of total energy with respect to the second 
state variable, the momentum p: 
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(P) (Parity of Total Energy with Respect to Momentum) For each F fixed, the 
total energy mapping is even in its second argument, namely, 

f(F,p) = f(F,-p). (3.18) 

Note that, due to (3.11), parity of total energy renders the velocity odd: 

v(F,p) = -v(F,-p). (3.19) 
The next two propositions state our main representation results. 

Proposition 3 (Constitutive representation of velocity) Let the constitutive 
mappings v and f satisfy, respectively, the variance and parity requirements (G) 
and (P). Then v has the representation 

v(F,p)=Vp, (3.20) 

where the symmetric and invertible tensor V depends at most on x. 
Proof. By assumption (G), 

9 F v(F,p + d) = 9 F v(F,p), (3.21) 
9 p v(F,p + d) = 9 p v(F,p), (3.22) 

so that, both 9pv and <9 p v do not depend on p. Relation (3.21) then yields the 
preliminary representation 

v(F,p)=V(F)p + v(F), (3.23) 

which, by consequence (3.18) of assumption (P), reduces to 

v(F,p)=V(F)p. (3.24) 

Consistency with (3.20) requires that V(F) = V, with V depending at most on x. 
Moreover, the normality assumption (N) guarantees that the tensor V is invertible. 
Finally, it follows from (3.11) and (3.23) that V is symmetric. □ 

Proposition 4 (Constitutive representation of energy) If Propositions (1)- 
(3) hold, then the total energy f(F,p) has the representation 

f(F,p)=/c(p)+a(F), £( P ) = ip-Vp, (3.25) 

where k(p) is interpreted as the kinetic energy and <r(F) as the stored energy. 

Proof. By (3.12), (3.13), and (3.19), we have that 

agr(F,p) = 0, (3.26) 
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which yields ( |3.25| )i. The representation ( |3.25 )2 then follows from (3.11) and again 
(3.19). □ 

Because of the presence of the tensor V, the expressions for v and k are slightly 
more general than the usual ones. As expected, the stored energy determines the 
stress, since ( 3.13 ) and ( 3.25| )i immediately yield 



S(F) = <9 f <t(F) 



(3.27) 



We note that, since (the velocity mapping is independent of F and) the stress 
mapping is independent of p, the ellipticity condition (E) can be phrased in terms 
of S(F,p). 



3.3. Constitutive representations and hyperbolicity 

We have proposed a set of reasoned assumptions under which the general constitu- 
tive relati ons Q2.18| ), ( |2.19D , and ( |2.20p take, respectively, the forms Q3.20| ), Q3.27| ), 
and d3~25|) . The latter group of relations implies hyperbolicity of the elasticity sys- 
tem, provided further assumptions are introduced. Such hyperbolicity — which, as 
is well-known, cannot be strict — is described precisely in our last proposition. 
As a premise, we note that, for 



§ = d F S = d^a 



the ellipticity matrix can be written as a linear transformation of vectors: 
E(w)u = (8[u <8> w])w, w,ueV, |w| = 1. 



(3.28) 



(3.29) 



Due to ( 3.28| ), strong ellipticity of the stress mapping S (that is, positivity of E(w) 
for each unit vector w) is equivalent to strict (rank 1 ) - convexity of the stored-energy 
mapping a. 

Proposition 5 (Hyperbolicity of Elasticity in Conservation Format) Let 

V = p i, where p > is interpreted as the referential mass density. More- 
over, let the matrix E(w) be positive for each unit vector w. Then, for each w, the 



(12x 12)-matrix that now corresponds to the general hyperbolicity matrix (2.3) has 
zero as proper value of geometric multiplicity six, and has, moreover, exactly three 
independent proper vectors corresponding to non-null proper numbers and giving 
the possible directions of wave propagation. 8 



Proof. Under the assumptions, for each fixed w, the problem we study can be 
given the form 

M (w) V = XV, (3.30) 



S The occurrence of zero as a proper number with high multiplicity is dwp to the involutions (2.9) 
which render the elasticity system not strictly hyperbolic (cf. DafermosQ). 
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where the entries of the (4 x 4)-blockmatrix M that satisfies ( |2.3| ) are the following 
tensors: 

0, 



» Q/3 =1*1144 
^a4 = -p _1 W (g) C Q , 
l4 Q U = — (§[u ® C Q ])w, 



{a, = 1,2, 3), 



u e V, 



(3.31) 



and where, given a tensor Z and a vector z, the entries of the 4-blockvector V are 
the vectors 

v Q = Zc Q (a = 1,2,3), V4 = z. (3.32) 
With ( 3.3l| ) the equation ( 3.30 ) yields the system 



M Q 4V4 = Av Q , 

3 

^M 4ct v Q = Av 4 , 



a=l 



or rather, more explicitly, 



— p 1 z ® w = AZ, 
-(S[Z])w = Az. 



(3.33) 
(3.34) 



(3.35) 
(3.36) 



For A = 0, (3.35) implies that z = 0; moreover, choosing C3 = w, we can 
represent Z as 

Z = vi ® ci + v 2 ® c 2 + v 3 ® w; (3.37) 
consequently, ( |3.36| ) becomes 



(S[v3 ® w])w = E(w)v3 = — (S[vi (8 Ci + V2 ® C2])w, 



(3.38) 



whence, in view of the assumed positivity of E(w), we conclude that Z belongs to 
a six-dimensional subspace of the space of all tensors. Moreover, assuming provi- 
sionally that A ^ 0, one finds by inserting (3.35) into ( 3.36) that 

E(w)z = /iz, p = p\ 2 . (3.39) 

Hence, the positivity of E(w) implies that there are positive proper numbers p and 
exactly three independent proper vectors. □ 



We remark that the involution condition (2.6) — which we do not use — implies that 
the proper space associated to the proper number A = is six-dimensional. 
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